Plasma resonance at low magnetic fields as a probe of vortex line meandering in 

layered superconductors 
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We consider the magnetic field dependence of the plasma resonance frequency in pristine and in 
irradiated Bi2Sr2CaCu2 08 crystals near T c . At low magnetic fields we relate linear in field corrections 
to the plasma frequency to the average distance between the pancake vortices in the neighboring 
layers (wandering length). We calculate the wandering length in the case of thermal wiggling of 
vortex lines, taking into account both Josephson and magnetic interlayer coupling of pancakes. 
Analyzing experimental data, we found that (i) the wandering length becomes comparable with the 
London penetration depth near T c and (ii) at small melting fields (< 20 G) the wandering length 
does not change much at the melting transition. This shows existence of the line liquid phase in this 
field range. We also found that pinning by columnar defects affects weakly the field dependence of 
the plasma resonance frequency near T c . 



Josephson plasma resonance (JPR) measurements 
in highly anisotropic layered superconductors provide 
unique information on the interlayer Josephson coupling 
and on the effect of pancake vortices on this coupling. 
The squared c-axis plasma resonance frequency, lo 2 , is 
proportional to the average interlayer Josephson energy, 



lo 2 cx Jq (cos (r)), where Jq is the Josephson critical 



current, </? n! „ +1 (r) is the gauge-invariant phase difference 
between layers n and n + 1 and r is the in-plane coordi- 
nate. Here (...) means average over thermal disorder and 
pinning. Thermal fluctuations and uncorrelated pinning 
lead to misalignment of pancake vortices induced by the 
magnetic field applied along the c axis. Misalignment re- 
sults in nonzero phase difference and in the suppression 
of Josephson coupling and plasma frequency!!! Thus, the 
uj p dependence on the c-axis magnetic field measures the 
c-axis correlations of pancakes in the vortex state. 

The JPR measurements performed in the liquid vor- 
tex phase at relatively high magnetic fields, B > Bj = 
<i>o/Aj revealed that the plasma frequency drops approx- 
imately as Here Xj = 7s is the Josephson length, 
7 is the anisotropy ratio and s is the interlayer dis- 
tance. The above dependence is characteristic Jar the 
pancake liquid weakly correlated along the c axis. Era Here 
a pancake in a given layer is shifted by a distance of 
the order of vortex spacing a = {Qq/B) 1 / 2 <C Xj from 
the nearest pancake in the neighboring layer. Thus at 
high fields many pancake vortices contribute to the sup- 
pression of the phase difference at a given point r, be- 
cause Xj determines the decay length for the phase dif- 
ference induced by misaligned pancakes of a given vor- 
tex linc.El In contrast, in the vortex solid a lattice of 
vortex lines .forms as shown by neutron scattering- and 
/i + SR data. El JPRj measurements in Bi2Sr2CaCu20s_5 
(Bi-2212) crystalsLI u have shown that in the fields above 
20 Oe the interlayer phase coherence changes drastically 
at the transition line [(cos (p„ jra +i(r)) jumps], implying 



the decoupling nature of the first-order melting transi- 
tion, see discussion in Ref. EoL 

In this Letter we focus on the low magnetic fields, 
B < Bj, near T c regime. Here the intervortex dis- 
tance is much larger than Aj, and the Josephson cou- 
pling in the region occupied by a given vortex is not sup- 
pressed by other vortices. In such a single vortex regime 
the Josephson energy increases linearly with displace- 
ments of nearest pancakes in neighboring layers, when 
|u JV „ + i| = |r n - r n+ i| > Xj, see Ref. | (here r„ is the 
coordinate of a pancake in the layer n). This leads to 
the confinement of pancakes in neighboring layers, and 
c-axis correlated pancakes (i.e. vortex lines) may be pre- 
served above the melting transition. We show that it is 
the linear decrease of oj 2 with B that characterizes such 
a vortex state. The linear dependence was observed ex- 
perimentally in Refs. in both solid and liquid vortex 
states in Bi-2212 crystals in fields below 20 Oe near T c 
providing evidence for a line structure of the vortex liquid 
state at low fields. 

We calculate the plasma frequency at low magnetic 
fields B, and near T c assuming that only vortices induced 
by the applied magnetic field B || c contribute to the field 
dependence of u> p . We, thus, ignore the contribution of 
thermally excited vortices and antivortices to the field 
dependence of the plasma frequency. 

The JER absorption is described by a simplified 



equatiorO for small oscillations of the phase difference 
ip' n n+1 (r,uj) induced by an external microwave electric 
field with the amplitude T> and the frequency u) applied 
along the c axis: 



uj{uj + iT c ) 



X 2 jLX7 2 r - V„(r) 



Vn.n+l 



iluT> 
4ttJ ' 



(1) 



Here V„(r) = cos(^ n! „ + i(r) — 1 is the effective poten- 
tial, (r) is the phase difference induced by vortices 
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misaligned due to thermal fluctuations and pinning in 
the absence of a microwave field. In Eq. ([j]) we neglect 
the time variations of </2 n)n +i (r, t) because the plasma 
frequency is much higher than the characteristic frequen- 
cies of vortex fluctuations, see below. Further, u>o(T) — 
c/^/eoX c (T) is the zero field plasma frequencjO, eo is the 
high frequency dielectric constant, X a b and A c = "f\ a b 
are the components of the London penetration depth, 
Ej = Eq/X'j is the Josephson energy per unit area and 
E = s^l/16n 3 X 2 b is the characteristic pancake energy. 
The inductive matrix L is defined as LA n — J2 m L n mA m 
with L nm ps (A a {,/2s) exp (— \n — m\s/X a b)- The parame- 
ter r c = 47rcT c /eo describes dissipation due to quasipar- 
ticles. Here a c is the c-axis quasiparticle conductivity in 
the superconducting state. Practically it coincides with 
conductivity right above T c . 

The absorption in the uniform AC electric field is de- 
fined by the imaginary part of the inverse dielectric func- 
tion 



Im 



— = -y 

e(u>) e j-f 



{u 2 - UJ 2 a ) 2 + UJ 2 T 2 C ' W 



where E a = 1 — w^/wq and ^ r an (r) are the eigenvalues 
and eigenfunctions of the operator — A 2 LV 2 + V„(r). 

Consider magnetic fields B <C $o I^X 2 ^, Bj (single 
vortex regime) . The phase difference near a given vortex 
line is induced by displacements of pancakes in neighbor- 
ing layers along this vortex line (see Fig. [j]). The poten- 
tial Vn(r) at distances u n , n +i <C r <C Xj is determined 
by the phase difference produced by nearest pancakes in 
neighboring layers n and n + 1 relatively displaced at 
distance u„ n +i: 



(r) « [r X Un, n+1 ]/j 



(3) 



At large distances r ^> Xj the potential drops exponen- 
tially and at small distances r < M n ,n+i it tends to a 
constant. This potential is attractive and there are lo- 
calized and delocalized states. At a > Xj,u n>n +i main 
contribution to absorption is coming from the most ho- 
mogeneous delocalized state. Such a state determines 
the center of JPR line uj p {B,T). Other states lead to 
inhomogeneous line broadening in addition to broaden- 
ing caused by quasiparticle dissipation described by T c . 
The latter mechanism of broadening dominates at low 
magnetic fields near T c . 

The strength of the potential with respect to the ki- 
netic term, AjLV 2 ., is characterized by the dimensionlcss 
parameter r^/X 2 which we assume to be small in the 
following calculations. Here r w is the elemental wander- 



ing length of vortex lines, r 2 



l n,n+lh 



Then we use 



perturbation theory with respect to the potential V„(r) 
to find the energy of the most uniform delocalized state. 
The unperturbed wave function of this state is given by 
a constant. The first order correction to the energy of 
this most homogeneous delocalized state is given by the 
space average of the potential, averaging is over r and 



n. This space average is equivalent to the thermal aver- 
age of cos (r) - 1 fa -if 2 n n+1 (r)/2. Using Eq. (g) 
we obtain a simple relation connecting field-induced sup- 
pression of the plasma frequency u p {B,T) with r w for 
the case r w < Xj: 

u 2 (T)-cj 2 JB,T) «„ +1 (r)) nBr. 



ln^. (4) 

2$o r w y ' 



The relation (Q) is very general and does not depend 
on the mechanism of the vortex wandering. It allows 
one to extract r 2 ^ from the plasma resonance measure- 
ments. The field dependence of the resonance tempera- 
ture, T r (B), is determined by the equation iv 2 (B,T r ) — 
lo 2 . According to Eq. (Q) this gives a linear dependence 
at small fields, T r (B) ps T r (0) + {dT r /dB)B, and r 2 w is 
directly related to the slope of this dependence 



Xj _ 2$ d\noj 2 (T) 



An— = 



dT 



dT r 
~dB 



(5) 



We now calculate when wandering of the vortex lines 
is caused by thermal fluctuations. In the single vortex 
regime is determined by the wandering energy con- 
sisting of the Josephson and magnetic contributions, 



E 



-n+l 



(6) 



where £\j P3 ($q/(47tA c ) 2 ) \n{Xj/r w ) is the line ten- 
sion due to the Josephson coupling and wm ~ 
{^ol {^X 2 ah ) 2 ) ln(X a b/r w ) is the effective cage potential, 
which appears due to strongly nonlocal magnetic mt j jr- . 
actions between pancake vortices in different layers. tM3 
Assuming Gaussian fluctuations we have 



' wT 



dq z 



4T(1 - cos q z s) 



/(C) = 



27r 2 (e\j / s 2 ) (1 — cosq 2 s) + wm 

c 



2sT 



1 + C + VTTC' 



/(C), 



(7) 



where the parameter £(T) = 4X 2 b (T)/X 2 describes the 
relative roles of the Josephson and magnetic interactions. 
Substituting this result into Eq. (0) we obtain 



u 2 (T)~u 2 (B,T) 



|-/(0 

Bo 



(8) 



with B = $g/167r 3 A 2 sT = Bj(E /T). We stress that 
this result of the single vortex regime is valid in both solid 
and liquid vortex states for B -C Bj, because in this field 
range wandering of lines at short scales does not change 
much at the melting point. The difference between these 
states appears only in the second order in the magnetic 
field. At B = the resonance occurs at the temperature 
T r (u). The slope of the curve B r (T) at small B r is 



dB r 



1 



/(Cr) 



dB 
dT 



4 = * (9) 



T=T r 



A 2 ; 
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The crossover region from the line liquid, <C a , to 
the pancake liquid, where ~ a 2 , takes place at the 
magnetic field B w min[7rSo/2/(C), -Bj]. 

To compare our calculation with experiment we plot 
in Fig. ^| dependence of dB r /dT vs reduced resonance 
temperature at B = 0, t r = (T c — T r )/T c , obtained in 
Ref. 0,| using different microwave frequencies (shown 
in the plot) for Bi-2212 with T c = 84.45j.jK. We also 
show data point obtained by Matsuda et alu for Bi-2212 
with close T c , T c = 85.7 K. To calculate the dependence 
dB r /dT from Eq. (||) we need dependencies A C (T) and 
^ab(T), which determine dependencies B (T) and C(^)- 
A C (T) was found from temperature dependence of luq at 
zero field, which we fit as uj {t)/27r m 033.5 GHz)t 32 , 
and taking e c = 11. Matsuda et aln obtained simi- 
lar temperature dependence of ojq using direct frequency 
scan. X a b(T) was obtained assuming temperature inde- 
pendent 7, which we adjusted to obtain the best agree- 
ment between the theoretical and experimental curves 
giving 7 = 480. Nonmonotonic temperature dependence 
of dB/dT r arises from competition between two factors 
in Eq. (pj): increase at low temperatures is due to the 
factor 1//(C) oc A~ b 2 at £ -C 1, and increase at temper- 
atures close to T c is due to nonlinearity of the depen- 
dence A~ 2 (T), which leads to the divergency of dBo/dT 
at T — ► T c , As one can see from the plot, our theory de- 
scribes satisfactorily the field dependence of oj p not very 
close to T c , i.e., in the region where the critical fluctua- 
tions are not very strong. The region of critical fluctua- 
tions is beyond applicability of our theory. 

We can now check validity of our approximations: the 
static approximation for the potential in Eq. ([!]) and 
the perturbation theory with respect to the potential. 
The maximum frequency of vortex fluctuations in the 
single vortex regime is wj; sa eu/s 2 i], where the vor- 
tex viscosity rj estimated from flux flow resistivity is in 
the interval 10 -6 — 10~ 7 g/cm-s, see Ref. 14 
s 2 A 



d> 2 



10" 

V167T 2 



< 0.1 near T c 



This gives 
Condition 



j/^j ~ 1 f° r applicability of the perturbation theory is 
satisfied at t > 0.02. 

Near T c The experimental curve T r (B) shows no 
change of slope when it crosses the melting line, H m 
1.3[Oe/K](T c - T)n This gives evidence that the param- 
eter r 2 does not change much at melting. Using Eq. (||) 
we estimate that, near melting line r^/a 2 < 0.1 at tem- 
peratures studied. This confirms the line structure of the 
vortex liquid above the melting line at low fields, though 
vortex lines wander over extended distances already in 
the solid phase due to high temperatures, see Fig. |]. The 
estimated value, r w « 1 /im at 77 K, is comparable with 
both Aj and A a h. In Bi-2212 crystals near optimal doping 
crossover to a pancake liquid occurs in the field interval 
~ 10 — 15 Oe. In less anisotropic high-T c materials one 
anticipates a larger region of the line liquid on the vortex 
phase diagram. 

Next we calculate the effect of columnar defects (CDs) 
on the parameter at high temperatures. Columnar 
defects always straighten vortex lines and reduce r w . At 



low temperatures this effect is very strong. Each vortex 
line is localized near one defect and its wandering is much 
smaller than in an unirradiated superconductor. At high 
temperatures lines start to distribute over a large number 
of defects and the effect of CDs progressively decreases. 
We consider the extreme case of very high temperatures 
when the effect of defects can be considered within per- 
turbation theory. This approach is applicable at tem- 
peratures higher than the pinning energy of pancakes by 
CD, T > ttE ln(b/£ a b), where b is CD radius and £ a f, 
is the superconducting correlation length. This situation 
corresponds to experiments. Qcl 
The free energy functional is 



(10) 



where U(r) = J2i ^( r — R-i) is the pinning potential of 
CDs, Ri are the positions of CDs, and V(r) = ttEq ln(l — 
b 2 /r 2 ) exp(=WA a b) is the pinning potential of the indi- 
vidual CD.E3 Expanding with respect to disorder up to 



second order terms we obtain the correction, (- 
the zero order term, r^ T , due to pinning by CDs: 



wD 



), to 



' wD 
S 



Tdei 



r — r 

w 

—Y. 



wT 



(11) 



J 



[(K(r m -ro)) -((K(r' -r )) )' 



where (. . .) stands for statistical average for a system 
without disorder and K(r) is the correlation function of 
disorder, 



K(v< -v) = (U(r')U(r)) D - (U(r)) 2 D . 



(12) 



When r w is much larger than the distance between 
columns we can transform Eq. (11) to a simpler form 



' wD 



2sK d 
T de u 



E 



1 



ro) 2 > 



l 0/0 



(13) 



where K = J dvK{v) = [it 2 b 2 E Q ln(X ab /b)] 2 . Here 
n<f, is the concentration of CDs. The correction is de- 
termined by the lateral line displacement ((r m — ro) 2 ) 
which we calculate as 



dq z 



4T(1 — cosmq z s) 



2ir 2 (eu/s 2 ) (1 - cos<j z s) + w M 

C+2-2VCTT 



<(r m -r ) 2 ) Q = 

_ 4sT v (1 - v m ) 
"717 (1-v 2 ) ' "~ C ' 

Combining Eqs. (^3|) and ( |l4| ) we finally obtain for 

2 2K 



(14) 



vD 



uD 



y2 



■fn(v). 



(15) 



Here the dimensionless function 



f D (v) 



= jl-v) d 
1 + v dv 



(l + vf- 1 ' 

I 1 - V ^ v~ n - 1 

\ n=l 
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has the limits /d(0) = 1 and Jd(v) ~ 2.16- 
v — > 1. Comparing this equation with Eq. 



■21n(l-u) at 
m) we obtain 



' wT 



In" 



A a £ 



T 



hi 



(16) 



In the crystals studied Aj f=s 1 /im, A ab (0) w 2000 A, 
b ~ 70 A. For the temperature range T > 77 K explored 
in Refs. |LB| by the plasma resonance this gives a very 
small correction (r 2 n/ r wT ~ 10~ 3 )- Therefore, the effect 
of pinning by CDs on the field dependence of the plasma 
frequency near T c is negligible. 

It was found in Refs. that after irradiation with the 
matching field = <&o"-0 ~ 1 T the value (dB/dT r ) B ^ 
near T c increases about two times in comparison with 
that in pristine crystals. As was estimated above, pinning 
due to CDs cannot give such a strong effect. One may 
think that irradiation reduces the value of the anisotropy 
parameter 7, probably due to damage of the crystal struc- 
ture around the heavy ion tracks. This assumption is 
consistent with recent measurements of the Josephs©*! 
current in irradiated Bi-2212 mesas by Yurgens et alx3 
It was found that irradiation approximately doubles the 
Josephson current at zero field. 

In conclusion, we have calculated the field dependence 
of the JPR frequency in the single vortex regime at 
low magnetic fields near T c and demonstrated that the 
JPR provides a direct probe for meandering of individ- 
ual lines. We have shown that the JPR data in highly 
anisotropic Bi-2212 crystals give evidence that at high 
magnetic fields B ^> Bj pancakes are uncorrelated along 
the c axis in the vortex liquid (pancake liquid), while at 
lower fields, B < Bj, pancakes form vortex lines (line 
liquid). These lines, however, strongly meander in both 
solid and liquid vortex states due to thermal fluctuations. 
We have shown also that JPR data provide evidence that 
irradiation by heavy ions causes a significant decrease of 
the effective anisotropy. 
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FIG. 1. Meandering of pancakes along the vortex line in 
the single vortex regime at low magnetic fields in highly 
anisotropic layered superconductors. 
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FIG. 2. Comparison of the experimental temperature de- 
pendence of dB r /dT obtained in Ref. |^ using different mi- 
crowave frequencies with theoretical dependence from Eq. (|j) 
(see text). We also show data point from Ref. |^ obtained for 
Bi-2212 with slightly different T c . 
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